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Abstract 

A kinetic equation for the joint probability distribution for fixed values of the 
classical action, momentum and density has been derived. Further the hydro- 
dynamic equations of continuity and balance of momentum density have been 
transformed into a Schroedinger-like equation, describing particle motion in an 
effective electro-magnetic field and an expression for the gauge potentials has 
been obtained. 



1 Introduction 



The dynamics of particles in accelerators and storage rings is usually studied on the basis of 
deterministic causal tools, provided by Hamilton's equations of motion. There are, however, 
a number of instances where such a description fails or may not be adequate. 

The beam circulating in an accelerator may be generally considered as an ensemble of 
nonlinear oscillators. Even in the case when the beam is not dominated by space charge 
these oscillators are weakly coupled at least linearly. This coupling is due to direct in- 
teraction between particles, thus revealing the discrete nature of their correlations and/or 
an interaction between them via the surroundings. In an experiment involving macroscopic 
measuring devices the observed quantities are a restricted number of variables characterizing 
the macroscopic state of the beam. Since the beam consists of a large number of randomly 
moving particles the macroscopic quantities are subject to deviations around certain mean 
values. These deviations appearing to an observer as random events are called fluctuations. 

Therefore the particle beam propagates in response to fluctuations, automatically implied 
by the existence of many degrees of freedom. Fluctuating contributions remain small in 
comparison with the macroscopic quantities for systems in thermodynamic equilibrium 
except at critical points. When a certain dynamic (or plasma) instability is encountered, 
fluctuations are expected to grow considerably. As a result the corresponding macroscopic 
evolution exhibits abrupt changes of various thermodynamic parameters. 

In the light of the above considerations we analyze in the present paper the motion of 
a test particle in the bath provided by the rest of the beam. Microscopically, each particle 
feels the fluctuating field (due to electro-magnetic and other possible interactions) produced 
by all the other particles in the beam, and therefore it constantly undergoes Brownian 
motion in phase space. The coefficients in the resulting Fokker-Planck equation contain the 
fluctuation spectrum of the interparticle interactions 0, We are not going to calculate 
these coefficients explicitly here (we hope to do so in a forthcoming paper), but rather we 
study the motion of a test particle suspended in a random inhomogeneous medium under 
the action of external forces. The statistical properties of the medium comprising the rest 
of the beam are characterized by a random velocity field (which may be regarded also as 
a fluctuating vector electro- magnetic potential) and a random potential fleld. Moreover, 
we consider the beam fluid inviscid. The latter restriction is not essential for we presume 
friction and other sources of dissipation (such as synchrotron radiation) could be without 
effort implemented in the development presented here. 

Recently a thermal wave model for relativistic charged particle beam propagation, build- 
ing on remarkable analogies between particle beam optics and non relativistic quantum 
mechanics has been proposed 0]. The conjectured in Reference 4 Schroedinger-like equa- 
tion for the transverse motion of azimuthally symmetric beams has been derived in the 
framework of Nelson's stochastic mechanics. Further development of the model suitable to 
cover the problem of asymmetric beam propagation in accelerators can be found in p. In 
the present paper we recover the Nelson's scheme of stochastic mechanics for particle beams 
from a different point of view. 
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2 The Model of Random Beam Propagation 



The classical motion of charged particles in an accelerator is described usually with respect 
to a comoving with the beam reference frame. Not taking into account chromatic effects 
the dynamics in the longitudinal direction can be decoupled from the dynamics in a plane 
transversal to the orbit. Then the evolution of the beam in 6D phase space is governed by 
the Hamiltonian [^: 

H{^,p;e) = ^+U{^;e), (2.1) 

where x = {xi,X2,X3), p = (^1,^25^3) and 9 is the azimuthal angle, commonly used in 
accelerator theory as an independent variable playing the role of time. The quantity R in 
equation (2.1) denotes the mean radius of the machine. The variables X3, ps constitute a 
canonical conjugate pair, responsible for the longitudinal motion of the beam 

X3 = ^sign{lC) ; ps = -r]J\}C\sign{}C), (2.2) 

l/CI 



where rj is the deviation of the actual energy E of the particle under consideration with 
respect to the energy Eg of the synchronous particle 

I E-E, ( Vs\ 

^ = ^ — F — ; \Ps = — ] (2.3) 



K E, 

and a is the displacement of the longitudinal position of the particle with respect to the 
longitudinal position of the synchronous particle 

a = s- (3sCt. (2.4) 
The quantity K, is the auto-phasing coefficient (phase slip factor), related to the momentum 



compaction factor au through the equation 



/C = «M-^ ; |7s= , ^ |. (2.5) 



Ts V ^1 - 

The beam propagation in the plane transversal to the particle orbit is described by the 
canonical conjugate pairs 

Xk = xu-7]E>k{e) ; Pk=Pk-^Dk{0) (A; = 1,2). (2.6) 

In equation (2.6) Xk is the actual position of our particle in the transversal plane and pk 
is the canonical conjugate momentum scaled by the total momentum ps = nioPsIsC of the 
synchronous particle. The function Dk{6) is the dispersion function defined as a solution of 
the equation 

Dk{0) + Gk{9)Dk{9) = R'Kk{9), (2.7) 
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where Gk{0) are the focusing strengths of the hnear machine in the two transverse directions, 
KkiO) is the local curvature of the orbit and the dot [as well as in equation (2.6)] stands for 
differentiation with respect to 9. 

The potential function W(x; 9) in equation (2.1) consists of two parts: 

U{^-9)=U,{x^,X2-9)+U,{x^-9), (2.8) 
where Ui,{xi,X2', 9) describes the transverse motion (betatron motion) and is given by 



Ub{xi,X2; 9) = — \Gi{9)xI + G2{9)xl\ + V(xi, xa; i 



(2.9) 



while Us{x^\ 9) is responsible for the longitudinal motion (synchrotron motion) and has the 
form: 



Us{x^]9) = -sign{IC) ^^-^cos 



uJ^J\IC\sign{}C) 



-X3 + ipo 



(2.10) 



In formula (2.10) AEq is the maximum energy gain per turn, u and ifo being the angular 
frequency and phase of the accelerating voltage respectively. 

From classical mechanics it is well-known that the Hamilton- Jacobi equation associated 
with the Hamiltonian (2.1) is 



dS{^;9) R 2, ,// 
— ^ + -p2(x; 9)+U{^; 



0, 



(2.11) 



where 



p(x;e)=V5(x;e). 



(2.12) 



For a given arbitrary integral of equation (2.11) a family of trajectories q(6') is generated 
that solve the first order (vector) differential equation 

cii9) = Rp[qi9y,9]. (2.13) 

Moreover the continuous distribution of trajectories with associated density ^(x; 9) obeys 
the continuity equation 



d0{x;9) 
39 



/?V-[^(x;^)p(x;^)] = 



(2.14) 



and in addition, taking the gradient of equation (2.11) we obtain an equation for p(x; 9) 

(9p 
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+ i?(p • V)p + VW(x; 9) = 0. 



(2.15) 



Thus the system (2.14), (2. 15) [or equivalently (2.11) and (2.14)] represents a closed set of 
equations, describing the Hamilton- Jacobi fluid as a mechanical system living in configura- 
tion space M. 
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Let us consider the motion of a test particle in the fluid comprised by the rest of the beam. 
No dissipative forces of Stokes type are present, as soon as we assume the beam fluid to be 
inviscid. However, the discrete nature of collisions between particles (intra-beam scattering) 
give rise to kinetic fluctuations in the one particle distribution function. As a consequence 
the gas-dynamic functions f)(x; 6) and p(x; 6) of the Hamilton- Jacobi fluid fluctuate as well. 
Fluctuating contributions to the one particle distribution function are generated also by the 
electro-magnetic interaction between particles in the beam. All this means that the beam 
is a real medium with a finite number of particles within a physically infinitesimal volume 
and substituting it by a continuous medium is not justified 0. External noise could be 
introduced into the beam from the surroundings (RF noise, fluctuations in the parameters 
of magnetic elements, etc.), which complexifies the physical picture additionally. 

Leaving more speculations aside we consider the motion of our test particle in a random 
inhomogeneous medium and random velocity field. The particle dynamics is governed by 
the Hamiltonian 

i7(x, p;e) = ^ + p- z{e) + W(x; e) + W(x; 6), (2.16) 
where z{6) is a random velocity field with formal correlation properties 

{z{e)) = ; {zk{e)zm{Oi)) = Rekm6{e - e,). (2.17) 

The quantity W(x; 6) is a random potential accounting for the fluctuation of the medium. 
In what follows we shall consider the random potential field W(x; 6) 5-correlated with zero 
mean value and correlation function A(x, Xi;^) 



[U{x;e)) = ; {U{x; e)U{xi; 61)) = A{x,xi; 6)5(6- Oi). (2.18) 

In the Hamiltonian (2.16) we have discarded the "constant" term z^{6)/2R for it does not 
give contribution to the dynamics. Moreover, e^m has the dimension of emittance and we 
call it thermal beam emittance tensor. 

We would like to note that the equation of random trajectories [replacing now equation 
(2.13)] 

q(0) = i?p[q(^);0]+z(0) (2.19) 
is in fact the equation for the characteristics of 

- = - + (fip + z).VC = (2.20) 

describing the mixing of concentrations C(x; 6) of different species in a random velocity field. 
Instead of equations (2.11), (2.14) and (2.15) we have now the system 

— + -p2 + p • zi6) + W(x; 6) + W(x; 6) = 0, (2.21a) 
^ + V-{[Rp + z{6M = 0, (2.216) 
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^ + [{Rp + z) ■ V]p + VW(x; 6) + VW(x; 0) = 0, (2.21c) 

Go 

which specifies the evolution law of the Hamilton- Jacobi fiuid with hydrodynamic Langevin 
sources. 

3 Kinetic Equation for the One- Point Probability Den- 
sity 

We define a joint probability density for fixed values of the classical action S'(x;6'), of the 
momentum p(x; 6) and the density of random trajectories ^(x; 6) as 

W{S,p,g\^;9) = {WriS,p, g\^; 9))^^^, (3.1) 

where denotes statistical average over the ensemble of realizations of the stochastic 

processes indicated. 

Note that now ^(x;^), p(x; 6') and f?(x; 0) are random functions [more precisely, func- 
tionals of the random velocity field z(^) and the random potential W(x; 9)] according to 
the system (2.21). A closed kinetic equation for the one-point probability density taking 
into account the gas-dynamic equations (2.21) can be derived by particular choice of the 
dependence of on the density f)(x; 9) [|ll|], that is 

Wr{S, p, f?|x; 9) = f?(x; 9)6[S{^- 9) - S]5[p(x; 9) - p\. (3.2) 

Differentiating equation (3.1) with respect to "time" 9 and using the gas-dynamic equations 
(2.21) it is straightforward to obtain the following kinetic equation 

dW /i?D^ \ dW 



-V- (ziy,) + — (wiy,) + Vp- (ly.vw). (3.3) 



d_ 

dS 

It is worthwhile to note that if we let Wr depend on f)(x; 9) through a generic function it 
will turn out that the only possibility to cancel terms proportional to V ■ p appearing in 
equation (3.3) is to allow be a linear function of ^(x;0). However, the kinetic equation 
for the one-point probability density (3.1) with an arbitrary dependence on the density of 
random trajectories can be found in a closed form if the Hessian matrix 

_ d'S{^, 9) 

rimn — a W'^) 



of the classical action is included in the joint probability density 0, [|T0| and the system 
(2.21) is supplemented with an equation for the quantity defined by equation (3.4). 

We still have not reached our final goal, since the right hand side of equation (3.3) contains 
the yet unknown correlators of the random velocity field, the random potential field and Wr- 
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In order to split the above mentioned correlations let us consider a generic functional Tl[F] 
of the random Gaussian tensor field F^^ ^ (r;6'). Then the following relation holds |ll|] - 

m 

(F,,...,,„(r)7^[F]) = / rf-ri(F,,....,„(r)F^„...,^„(rO)( ,p ^'^''^\ A (3.5) 

which is known as the Furutsu-Novikov formula. In (3.5) r collects all the continuous ar- 
guments of the random tensor field, 5/(5Fmi,...,m„(i") denotes the functional derivative with 
respect to the random field and summation over repeated indices is implied. To apply the 
Furutsu- Novikov formula (3.5) we need the functional derivatives of 5'(x;6'), ^(x; 6') and 
p(x; 6*) with respect to the random velocity field z(^) and the random potential W(x;6'). 
From equations (2.21a-c) it is easy to find 



——— = -Pki^-e) ; -J- - = -5x-Xi ; ~ 77 = 0, 3.6a 

0Zk{6) (5W(xi;6') (5VifeW(xi;6') 



6zk{e) dxk ' 5U{^i-e) ' 5VifcW(x 



(3.66) 



(5p„(x;^) dpm{x;9) 5p^{^]e) Spm{yi;0) \ to a \ 

— ; 777 — = 7 ; =U : = — OfcmOiX — Xi). (3.6c) 

Szkio) dxk 5U{^i-e) 5VifcW(xi;^) ^ I) y ) 

By virtue of (3.6), (2.17) and (2.18) we cast equation (3.3) into the form: 

= VfcV^iy 2~~dS^ — 2 — ^P^^P™^' ^^-^^^ 

where we have taken into account the expansion of the correlation function (2.18): 

A(x, y; 6) = A{9) + Bkie){xk - Vk) + ]fkra{0){xk - yk){x^m -ym) + ---. (3.8) 

Without loss of generality the first term in the Taylor expansion (3.8) of the correlation 
function can be taken equal to zero, since it does not contribute to the dynamics (it embeds 
the gauge properties of the random potential field and can be scaled to zero). Thus we finally 
arrive at the desired kinetic equation for the one- point probability density: 

= ^VuV^W + ^V,,V,^l^, (3.7) 
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The kinetic equation (3.7) is rather comphcated to be solved directly, so approximate meth- 
ods to analyze it should be involved. For that purpose let us integrate (3.7) over S, that is 
exclude the classical action from consideration. We get 



^""^ I D V7 V77/ V7 R^km^ ^ , Ckm{()) 

— + Rp -Vw -VU ■ VpW = -^VkVmW + — - — Vp^Vp^w, 

where 

w{p,g\x; 9) = J dSW{S,p, g\x; 9). 
If we further integrate equation (3.9) over p we obtain 



39 



+ V- (f?)v 



'{+) 



Re 



km 



VfcV™(^) = 0, 



where 



(^(x; 9)) = J dpw{p, ^|x; 9), 

{g{x; 0))v(+)(x; = R j (ippw(p,^|x; 9] 
Defining the osmotic velocity u(x; 9) according to the Pick's law 

(^(x;^))«,(x;^) = -^V^(^(x;^)) 
one can write (3.11) in the form of a continuity equation 
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+ V-((^?)v) = 0, 



where 



v(x; 0) = V(+)(x; ^) + u(x; 9, 
is the current velocity. Next we introduce the stress tensor 



nmn(x; 9) = j dpp.mPnw{p, ^|x; i 



which consists of two parts 



The second term in equation (3.17) 

Gmni^; 9) = I dp\Rpm - f \RPn " 



U7(p, ^|x; 



(3.9) 



(3.10) 



(3.11) 



(3.12a) 



(3.126) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



is called the internal stress tensor. Multiplying the kinetic equation (3.9) by -Rp and inte- 
grating over p we obtain the transport equation for the momentum density 
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+ Va 



Re 



km 



-R{g)VnK-Vkgkn, (3.19) 



or in alternative form 



de 



+ 



■ V 



Re 



km 



VfcVm'y(+)^ 



-RV nU — T^'^kGkn, 



(3.20a) 



where use has been made of equations (3.16-18) and (3.11), and the backward velocity field 

V(_)(x;0) = v(x;0) + u(x;^) (3.21) 

has been introduced. One can immediately recognize in the left hand side of equation (3.20a) 
the mean backward derivative 0, 0, |T^ of the forward velocity 

1 



I?(_)t;(+)„(x; 6) = -RVnU - J^^kOkn- (3.20) 
Perform now "time" inversion transformation in equation (3.20) according to the relations 



e-^e' = -e x(^) ^ x'(^') = x(e) ; v{e) -^V{e') = -v{e). (3.22) 

As a consequence of (3.22) one has 

QQ, = -^ ' V,. = V, ; = -V,. (3.23a) 

In addition the forward and backward velocities and mean derivatives are transformed as 
follows i 

V(±)(x; e) v'(±)(x'; 6') = -V(^)(x; 9), (3.236) 

^^(±)/'(x';^') = -^^(T)/(x;^)> (3-23c) 

where /(x;^) is a generic function. Since the internal stress tensor Q^n is a dynamic char- 
acteristic of motion under time inversion its divergence changes sign. This also follows from 
the particular form of the "collision integral" [the right hand side of the kinetic equation 
(3.7)]. Therefore from (3.20) with (3.22) and (3.23) in hand we obtain 



P(+)t;(_)„(x; 9) = -RVM + —V^Gkn- (3.24) 



1 

(7) 

Equations (3.20) and (3.24) provide us two opportunities. First, summing them up we 
express the transport equation for the momentum density in terms of the current and osmotic 
velocities as 
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_ + (v • V)v = ~RVU + (u ■ V)u - ^VfcV^u. (3.25) 

This is nothing else but the Nelson's stochastic generalization of Newton's law |T1|, [T^. 
Secondly, subtracting equations (3.20) and (3.24) we obtain an equation for the internal 
stress tensor to be determined, that is: 

+ (V • V)n„ - (U ■ V)f„ + ^^Y^VkVmVn = J^^kGkn- (3.26) 

In the isotropic case ekm = ^^km (see the next section) by virtue of the equation 

^ + Vi(ui.vi) = ^Vi(Vi-vO (3.27) 

relating the current and osmotic velocity (which is a direct consequence of Pick's law and 
the continuity equation) we arrive at the following expression for the internal stress tensor: 

Gkl = -^{Q){^lkVln + VlnVlk), (3.28) 

where [compare with equation (4.2a)] 

= (MgM'^) = MkmMlnQ^n. (3.29) 



' kl 

Transforming back (3.28) to the original coordinates we obtain 

Gkn — 2 i^km^ m'^n "I" ^nm^ k'^m)- (3.30) 

Resuming the results of the present section it should be mentioned that the continuity 
equation (3.11) and the transport equation for the momentum density (3.19) are equivalent 
to Nelson's scheme of stochastic mechanics. Let us also note that the characteristics of 
the fluctuating beam medium, embedded in the random potential (2.18) do not enter the 
simple hydrodynamic approximation procedure adopted here up to the second moment. It 
remains, however to analyze the corrections to the evolution law of the Madelung fluid by 
taking into account the balance equation for the kinetic energy density. This can be done 
by employing more complete closure techniques to accomplish the transition between kinetic 
and hydrodynamic description. 



4 The Schroedinger-Like Equation 

Our starting point is the system of partial differential equations (3.13), (3.14) and (3.25) 
derived in the preceding section, which in fact represents the set of equations describing the 
evolution of the Madelung fluid in stochastic mechanics ||, |1^. Following [|| we perform a 
coordinate transformation 

Xi = My. {Xin = MnmXm), (4.1) 
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such that the transformed emittance tensor 



is proportional to the unit tensor 6ki 



kl 



mn 



(4.2a) 



e'ki = eSki (4.26) 

by a factor e, where {■ ■ ■)'^ denotes matrix transposition. The scaling factor can be chosen 
any of the eigenvalues {k = 1, 2, 3) of the original emittance tensor en- Provided eki is 
symmetric, the matrix A4 has the following structure 



M = AO, 

where O is an orthogonal matrix and A is the diagonal matrix 



(4.3) 



A, 



kl 



Jkl- 



Furthermore, the transformed current and osmotic velocities are |T6| 

vi = Mv ; ui = Mu, 
while the probability density in the new random coordinates is 

^1 



detM 



Then the transformed equations of stochastic mechanics read as 



ae 



r + Vi-(f?ivi) = o, 



QiUi = — Y^iQu 



dvi 



Re. 



where 



(vi • Vi)vi = -ReVM + (ui ■ Vi)ui - y^i^i' 



1 



(V,)„ = -Vi„. 

We are looking now for a Schroedinger-like equation of the type 
equivalent to the system (4.6) through the well-known de Broglie ansatz 



V^(xi;^) = \/^i(xi;^) 



exp 



(4.4) 

(4.5a) 
(4.56) 

(4.6a) 
(4.66) 
(4.6c) 

(4.7) 

(4.8) 

(4.9) 
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where H is a second order differential operator with known (constant) coefficients in front 
of the second derivatives. The basic requirement we impose on the operator H is to be 
Hermitian 



(4.10) 

which defines it (as can be easily shown) up to a generic scalar and vector functions. Without 
loss of generality one can write 



H = ^[iReVi + A(xi; 9)]' + <l>(xi; 6), 



(4.11) 



where the vector function A(xi;0) and the scalar function $(xi;^) define some effective 
electro-magnetic field. Substitution of the ansatz (4.9) into equation (4.8) followed by sepa- 
ration of terms by real and imaginary part yields: 



vi = Vi5 - A, 



(4.12a) 



^ + (vi ■ Vi)vi = E + vi X B + (ui • Vi)ui - 



where [17 



E 



_Vi$ 



B = Vi X A. 



(4.126) 



(4.13) 



Comparing equation (4.12b) with equation (4.6c) we conclude that the transformed external 
force —ReVJA equals the force produced by the effective electro-magnetic field: 



-ReVJA = E + vi X B. 



(4.14) 



In order to find the electro-magnetic potentials A(xi;^^) and $(xi;^^) we note that the 
Schroedinger-like equation (4.8) is gauge invariant under local phase transformation 



V'i(xi;6l) = ^(xi;6l)exp 



Re 



>Si(xi; I 



(4.15) 



This implies that the equation for the new wave function ?/'i(xi; 6) has the same structure 
as (4.8) with 



A(xi; 



Ai(xi; I 



A(xi;0) + Vi5i(xi; 



(4.16a) 



$(xi;0)— >$i(xi;0) = $(xi; 
Moreover, equation (4.14) written in the form 

dA 



dSi (xi 



36 



86 



+ (Vi5 - A) X (Vi X A) 



(4.166) 



(4.17) 
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is gauge invariant under (4.15) with 

5(xi; 6) 5'(xi; 6) = 5(xi; 6) + 5i(xi; 6). (4.18) 
Choosing iSi(xi; 9) = — 5(xi; 9) we obtain the Euler equation 

-ReVM + Vi'I'i = - Ai X (Vi X Ai) (4.19) 

for the gauge electro-magnetic potentials Ai(xi; 9) and $i(xi; 9). 

According to (2.8) - (2.10) the external potential W(x; 9) entering the Hamiltonian (2.1) 
can be specified as 

W(x;0) =Wo(x;0) + V(x;^). (4.20) 
The term Wo(x; 9) governs the linear motion and can be written in the form 

ZYo(x;0) = Vg(^)x, (4.21) 

where the matrix G(^) is symmetric in general, while V(x; 9) is a sum of all nonlinear terms. 
Further, we split the electric potential $i(xi; 9) into two parts according to the relation 

$i(xi; 9) = <l>o(xi; 9) + <l>;(xi; 9), (4.22) 

where 

$o(xi; 9) = |xf Gi(^)xi ; Gi(^) = MG{9)M-\ (4.23a) 

<|.;(xi;e) = -iA?(xi;0). (4.236) 
Equation (4.19) takes now the form 

-ReV^V = + (Ai ■ Vi)Ai, (4.24) 

which transformed back to the original coordinates x reads as 

(9 A' 

-i?VV = + (A' • V) A' (a' = MAi) . (4.25) 

o9 ^ ^ 

Equation (4.12a) suggests an alternative interpretation of the vector potential A(xi; 9) as the 
vortex part of the current velocity Vi(xi; 9). Taking into account (4.16a) and the particular 
choice of the gauge phase iSi(xi;^) one can expect that A'(x; ^) will be vortex-free in the 
original coordinates 

A'(x; 9) = -i?V<^(x; 9), (4.26) 
where ip(yi; 9) is the velocity potential |jl8|. The first integral of the equation (4.25) is then 
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^ + f|V^(x;0)f + V(x;0)=jW. 



(4.27) 



Without loss of generality the generic function g{6) may be set equal to zero as a result of 
the uncertainty in the definition of the velocity potential (4.26). Equation (4.27) is noting 
else but the Hamilton-Jacobi equation (2.11) for the "classical action" (^(x;^), associated 
with the nonlinear part V(x; 9) of the external potential ZY(x; 9) only. 

Performing a second [similar to (4.15)] phase transformation according to 



V^2(xi;6') = V'i(xi;6')exp 



—Rip{^,;9) 



(4.28) 



wc obtain the gauge potentials A2(xi;^) and $2(xi;^) entering the Schroedinger equation 
for the wave function ip2{^i', 9). They are: 



A2(xi;^) = i?(l-^2)Vi<^(A^"^Xi;^), 
$2(xi; 9) = $o(xi; 9) + RV{M-W. o) + 



(4.29a) 



The Schroedinger-like equation called upon to replace equation (4.8) reads as 
dM^i;9) 



iRe- 
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^[2i?eVi + A2(xi;^^)]V2(xi; 



$2(xi;^)^2(xi;i 



Retracing the sequence of phase transformations (4.15) and (4.28) we find 

- Rip[M-^^i]9) 

A2(xi; 9) = A(xi; 9) - Vi [5(xi; 9) - Rif[M-W, 9) 



^(xi; 9) = ^2(xi; ^)exp \ — [5(xi; 



$2(xi; 9) = $(xi; 9) + ^ [5(xi; 9) - Rip[M-W, 9) 
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(4.2%) 
(4.30) 

(4.31a) 
(4.316) 

(4.31c) 



The relations (4.31) indicate the equivalence of the Schroedinger-like equations (4.8) and 
(4.30) up to a global phase transformation, defined by the constant in the coordinates and 
time phase 

C = 5(xi; 9) - Rifi^M'W, O) = const. (4.32) 

The anisotropy of the random velocity field (2.17) reflects on the appearance of the gauge 
electro-magnetic potentials. There are two cases in which the vector potential A2(xi;^) 
vanishes and the scalar potential $2(xi;^) is equal (up to a non essential factor R) to the 
external potential W(xi;^). The first case is when the external potential is the harmonic 
oscillator potential (V = 0) with generally time-dependent frequency, while the second is the 
isotropic case {ekm = ^hm)- 
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5 Concluding Remarks 



In the present paper we have studied the motion of a test particle in a random inhomogeneous 
medium comprised by the rest of the beam. As a result of the investigation performed we 
have shown that Nelson's scheme of stochastic mechanics for particle beams in the case of 
zero friction, is equivalent to hydrodynamic approximation in the kinetic equation for the one 
particle distribution function up to the second moment. Further, it has been pointed out that 
the hydrodynamic equations of continuity and momentum density can be transformed by a 
change of coordinates and dependent variables into a Schroedinger-like equation. Regardless 
of the type of the external forces one need to introduce a gauge electro-magnetic field. If the 
beam constitutes an isotropic medium (holding in the case of symmetric beams) the gauge 
vector potential vanishes and as a consequence the scalar potential is equal to the potential 
that accounts for the external force. 

The gauge transformation (4.16) is the well-known transformation in classical electro- 
magnetic theory []T^ introduced by Weyl, indicating a transition to alternative electro- 
magnetic potentials, which sometimes are easier to find compared to the original ones. Be- 
sides that, the transformed potentials define the same electro- magnetic field tensor. Taking 
into account this fact we have found the gauge electro-magnetic potentials explicitly, depend- 
ing on the solution of a Hamilton- Jacobi equation for the classical motion of the particle in 
the anharmonic part of the external potential. 

The beam circulating in an accelerator consists of a large number of particles. Obviously, 
all of them cannot be in the same micro-state. As a result the beam itself generates noise, 
which plays a role similar to the role of perturbation in stability theory. The essential 
difference is that here the perturbation is produced by the macroscopic system itself (in 
addition to the noise introduced from the surroundings). In the present work we have 
adopted a phenomenological approach to describe beam fluctuations. In this connection it 
remains to compute the statistical properties of the beam medium in terms of the fluctuation 
spectrum, which we hope to perform in a forthcoming paper. 
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